Nonlinear Effects in Multi-Photon Polaritonics 



A. A. Pervishko,! T.C.H. Liew/-^ V. M. Kovalev,^ I. G. Savenko,!'"^ and I. A. Shelykhi'" 

^Division of Physics and Applied Physics, Nanyang Technological University 637371, Singapore 
^Mediterranean Institute of Fundamental Physics, 31 via Appia Nuova, OOO4O Rome, Italy 
'Russian Acad Sci, Inst Semicond Phys, Siberian Branch, Novosibirsk 630090, Russia 
'^Science Institute, University of Iceland, Dunhagi-3, IS-107, Reykjavik, Iceland 

(Dated: February 14, 2013) 

We consider theoretically nonlinear effects in a semiconductor quantum well embedded inside a 
photonic microcavity. Two-photon absorption by a 2p exciton state is considered and investigated; 
the matrix element of two-photon absorption is calculated. We compute the emission spectrum 
of the sample and demonstrate that under coherent pumping the nonlinearity of the two photon 
absorption process gives rise to bistability. 
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Introduction. — Polaritonics is a rapidly developing 
branch of science lying at the intersection of the physics 
of semiconductors, quantum mechanics and nonlinear op- 
tics. The interest in this field has been stimulated by 
a series of fascinating experimental and theoretical dis- 
coveries, which demonstrate that a variety of quantum 
collective phenomena can be observed in polariton sys- 
tems including: Bose-Einstein condensation (BEG) [Il- 
ls] and superfluidity [4H6]; the Josephson effect [7]; and 
the formation of various topological excitations such as 
vortices [BHTO] and solitons [TTHT^ . as well as other spa- 
tial patterns [TSHI7]. While most fundamental studies 
have been performed at low temperature, experiments 
in GaN [18l [19] , ZnO [20] and organic [21] microcavities 
demonstrate devices working at room temperature. Now 
it becomes obvious, that polaritonics has a huge potential 
for technological applications [Hj. 

Polaritonics stems from the proposal that quantum mi- 
crocavities can be used as efficient sources of THz radi- 
ation [23j . which needs efficient pumping of the 2p ex- 
citonic state. It can not be excited under single pho- 
ton absorption due to optical selection rules. In one 
of the recent publications [IS], some theoretical aspects 
of THz radiation emission under two-photon excitation 
were considered. The consideration was based on the 
master equation approach. The two-photon absorption 
and THz emission were investigated as incoherent pro- 
cesses. This kind of consideration, however, is not always 
applicable. One can imagine the situation when the two- 
photon transition is tuned in resonance with the energy 
of the confined electromagnetic mode of a cavity. In this 
case, the coupling is of resonant character and the the- 
oretical description should be performed along another 
route. Theoretical consideration of this regime is the 
subject of the present paper. 

We consider a cavity which is subject to two-photon 
resonant excitation, in which takes place, pumping of 
the dark (2p) exciton state. Instead of considering the 
strong coupling between near resonant photons and Is 
excitons [2i3j, we uncover a different kind of exciton- 



polariton generated by the non-linear coupling of pho- 
ton pairs with 2p excitons. Similiar to the conventional 
case, a characteristic splitting of the photoluminescence 
spectrum is observed. However, the value of the Rabi 
splitting now becomes dependent on the intensity of the 
pump due to the nonlinear nature of the two photon ab- 
sorption. This nonlinearity also results in the onset of 
bistable behavior. 

Photonics with a 2p state and the manifold. — We con- 
sider the following geometry (Fig. 1): 





FIG. 1: Geometry of the structure. We consider a micro- 
cavity, which is made by Bragg mirrors, with the quantum 
well inside, where a 2p excitonic state can be excited by two 
photons with energy hhJ2c each. 

a quantum well (QW) with 2p excitonic transition hav- 
ing energy hujp placed inside a planar microcavity sup- 
porting the photonic mode of energy u!2c ~ 'j-'p/2. In 
this configuration the cavity mode can resonantly inter- 
act with the 2p excitonic state. The Hamiltonian of the 
system reads: 

7i = hujpp^p + huj 20.0^ a + g {pa) a) + p^aa) (1) 
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where the operators p, a correspond to 2p dark excitons 
and photons respectively and satisfy bosonic commuta- 
tion relations. The first two terms in Eq. ([T]) describe 
free excitons and photons and the last term corresponds 
to the resonant interaction between them, where g is the 
two photon-exciton coupling constant calculated in the 
'Supplementary Materials' appendix using second-order 
perturbation theory: 
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where S" is a quantization area of the sample; mo, m* 
and /I are the free electron, effective and reduced exciton 
masses, respectively; q is the elementary charge of the 
electron; Eg is the band gap of the material; w is the 
energy of the photon; _E„ are the eigenvalues of the 2D 
Hydrogen atom; Rn{r) and normalized radial 

and angular eigenfunctions of the 2D Hydrogen atom, 
respectively, and A is the vector potential, which can be 
writen in the dipole approximation: 



A ^eAo 
Ao^ 



h 

2eeowiS'' 



where eo and e are vacuum and material permitivities; 
L is the length of the cavity. 

Note, that the Hamiltonian ([I]) commutes with the ex- 
citations' number operator N = a^a + 2p^p, that means 
that the eigenvectors of these two operators coincide and 
the Hilbert space of the problem can be represented as 
a direct product of the manifolds corresponding to given 
numbers of TV. Each manifold consists of [-/V/2J basis 
vectors: |A^,0),|A^ - 2, 1), |0, 7V/2) for even TV and 
|iV,0),|TV - 2,1),....,|1,/V - 1/2) for odd N. Here the 
first number in kets corresponds to the number of pho- 
tons, and the second number to the number of 2p exci- 
tons. The matrix element of Hamiltonian ([T]) between 
any two vectors belonging to manifolds with different N 
is zero. This means that determination of the spectrum 
of the problem consists in diagonalization oi N x N ma- 
trices, each of which has a finite size. 

For = and A'' = 1 the result is trivial and the cor- 
responding energies are Eq — and Ei — hi02c respec- 
tively. The cases of TV = 2, 3 are also straightforward to 
consider, as the size of the reduced Hamiltonian for the 



manifold is 2x2 in these cases: 
tt 
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and the eigenenergies in the case of resonance 2u}2c = '^p 
are 

E2 = 2hujc ± gV2, E^ = "ihuic ± <?V6. 

Consideration of the states with larger numbers of pho- 
tons requires diagonalization of matrices of higher size. 
The energy spectrum is represented in Fig. 2. 
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FIG. 2: Energy levels of the QW with 2p excitonic transition 
placed inside the microcavity. We consider the case when the 
exciton energy is hujp and the photon energy is lj2c ~ ijjp/2. 
Possible transitions for — 1,2,3,4 manifolds for the ideal 
cavity are represented in the right-hand figure. 

Spectrum of emission. — If the mirrors of the cavity are 
not perfect and some photons tunnel through them, the 
system can emit in the outside world. The spectrum of 
emission is determined by the energy distances between 
the levels corresponding to N and A'^ — 1 manifolds and 
the intensities of the transitions will be proportional to 
squares of the matrix elements of the photon annihilation 
operator, /,/ - |(*°"*|a|*™)pp(*™), where 1^-™) is an 
initial wavefunction in the manifold with N cavity pho- 
tons and I^P""*) is a final wavefunction in the manifold 
with N-1 cavity photons. represents the proba- 

bility of occupation of the initial state |27]. Assuming a 
laser excitation of the system, this is given by the Poisson 
distribution [3S]. Fig. 3, shows the total photolumines- 
cence emission spectrum. 
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FIG. 3: Emission spectrum of the system. The intensity of 
transitions from individual levels is presented by the black 
lines. We consider the case when uj2c ~ 'j-'p/S — 0.71eV,g = 
0.025me1/ and the average number of the excitations in the 
system (A^) = 30 and the statistics of photons is Poissonian. 
Including a Lorentzian broadening of each transition line gives 
the double peak spectrum of emission shown by the red curve. 
One can define the value of the Rabi splitting in the system 
as the distance between the two peaks. 
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Then, one should take averages of the above equations 
thus passing from the operators to the mean values as 
p ^ (p) — Tr{pp}. We present a mean product of several 
operators by products of the mean values of the operators 
(e.g. {pd'') — (p)(a)*,(a^) — )■ (a)^). This approximation 
is relevant for coherent statistics. Finally, we phenomeno- 
logically introduce the lifetimes of the modes and exter- 
nal coherent pumping of cavity photons (by a coherent 
laser beam) with frequency uj. As a result, Eqs. (2,3) 
transform into the following system of nonlinear differ- 
ential equations: 
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where 7p^a = l/(2rp.a) with Tp.a being lifetimes of the 
modes, P is the amplitude of the coherent pump and 
f2p — g/h. Making a substitution d <— ae~"^*,p ^ 
p^~2iLjt^ in the stationary regime the equations read: 



It is well-known, that in the case of one photon and 
a resonant Is bright exciton in a standard semiconduc- 
tor microcavity the energy splitting does not depend on 
the occupation number and equals AfJ^^s — 2gs, if gg 
is the bright exciton- single photon coupling strength. 
Considering our system with two-photon absorption for 
large N {N >> 1), and substituting a, \/N one im- 
mediately obtains Aflu w 2g\/N . Accounting for finite 
lifetime of the photons and excitons, one obtains (29j : 
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Varying the number of the excitations, which can be 
done by increasing of the pump power, one can thus 
change the Rabi splitting in the system. Note that cou- 
pling constant g scales with the quantization area as 
g ^ (S*)^^/^, so g'^N depends on the concentration of 
the photons Ua = N/S only. 

Classical field equations. — To investigate the dynam- 
ics of the system under coherent pump with account for 
finite lifetimes of the photons and dark excitons in the 
cavity, the approach should be modified. We will use a 
system of coupled nonlinear equations for the dynamics of 
the macroscopic wavefunctions describing cavity photons 
and dark excitons. The latter can be obtained in the fol- 
lowing way. First, let us write the Heisenberg equations 
of motion for the operators p and d: 



{ihAp + "fp)p + iMlpO^ = 0, 
{ihAa + 7a)a + 2iMlppa* = ^iP, 

where Ap = L0p — 2uj, Aa = uja —uj, which can be reduced 
to: 



{ihAa + ^a)a + - 
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Wa + iP = 0, 



This result can also be written in terms of the real func- 
tions describing the number of the photons, Na = |ap 
and excitons, Np — |pp: 
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Np = 



CiNa+ C2NI 

g'Nl 



(7) 



where 
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Equation (7) is a cubic equation for Na- The results are 
represented in Fig. 4 for the following parameters corre- 
sponding to the system: 7a = 0.05 meV, jp = 0.01 meV, 
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FIG. 4: Considering the finite Hfetime of the photons and ex- 
ternal coherent pumping of the cavity, we observed a bistable 
behaviour of the number of excitons on intensity of the coher- 
ent pump. The parameters are varied for detuning of the sys- 
tem: (dashed line) ftA^ — 0.375meV,hAp = 0.75meV; (dot 
line) hAa = 0.25meV,ftAp = 0.5meV; (bold line) HAa = 
0.125meV,hAp = 0.25meV. The influence of the exciton- 
exciton interaction on the hysteresis curve is illustrated in the 
inset of the figure, where the cases of Qp=0 (dashed curve) 
and Op — 6EBa%/S (solid curve) are presented. The detun- 
ing parameters are hAa ~ 0.375meV, hAp — 0.75meV. 

and three different values of detunings and Ap. These 
curves reveal an S-shapcd dependence of the mode popu- 
lations on the intensity of the pump, which characterizes 
the phenomenon of bistabihty, which is well-known in 
semiconductor microcavities [30 36 . If the value of the 
pumping intensity lies in the bistable zone, the system 
can, in principle, occupy more than one possible state, 
and the particular choice of this state depends on the 
history of the evolution. 

Bistability is a fundamental ingredient of several de- 
vices based on semiconductor microcavities, such as 
memory elements [37], spin switches [3S], spin pat- 
terns [3ni HO]) and optical circuit designs [HJ 
However, there is substantial difference between the pre- 
vious works and the paradigm developed in our Letter. 
Indeed, in most of the previous approaches bistability re- 
sulted from the exciton-exciton (Coulomb and exchange) 
interaction, even in THz emitting systems [33]. Here, 
the bistability arises due to the nonlinearity of the two- 
photon absorption process. It should be noted that, 
still, excitons are present in the system and their inter- 
action can, in principle, change the physical behavior. 
The exciton-exciton interaction can be accounted for by 
means of the term: 

nint = -jP P PP^ 
which should be added to Hamiltonian (1). Then, the 



equations of motion change to: 

dp , -0-2 « 

at n n 

and in the stationary regime we obtain the following so- 
lution: 

4((7a7p - h'^AaAp)np - hAaapU^p) 

To define the constant ap, we use the formula Up « 
GEbci^/ S, where Eb is the exciton binding energy, qb is 
the exciton Bohr radius and 5* is the laser spot area [44] . 
In QWs based on GaAs-based alloys, we can use the pa- 
rameters: Eb ~ 19.2 meV and as — 5.8 nm. The result 
is presented in Fig. 4. 

As one can see the exciton interaction term changes 
the shape of the bistable behavior of the system. It be- 
comes more extended on intensity of the coherent pump. 
The density of the excitons compared with the non- 
interacting case decreases as well. We therefore conclude 
that the influence of the nonlinearity of the two-photon 
absorption on the bistability of the system is comparable 
with the contribution of exciton-exciton interactions. 

Conclusions. — We developed a theoretical approach 
for the description of strong light-matter interaction in 
a system with two-photon absorption processes. We cal- 
culated the photoluminescence spectrum, where a split- 
ting occurs due to the coupling between photons and 2p 
dark-excitons in the two photon absorption process. We 
developed the mean-field equations, valid under coherent 
excitation of the system, and demonstrated that the non- 
linearity of the two-photon absorption process gives rise 
to the phenomenon of bistability. The two-photon to 2p 
exciton matrix element was accurately calculated. This 
formalism is not only directly applicable for device con- 
struction, but also highlights the nonlinear nature of the 
system, which may be useful for the realization of quan- 
tum optical devices under coherent pumping [JS] US] . 
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